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ABSTRACT 


This is a preliminary investigation, both 
theoretical and experimental on the wing-body interference 
of a class of bodies, wings and their •ombination. The 
formulation of the problem holds good both in supersonic 
and subsonic flov/ regimes based on linearised equations. 

A case of a wing alone and body alone are studied both 
experimentally and theoretically in a subsonic flow regime. 
A different wing-body combination has been experimentally 
investigated for the combined aerodynamic characteristics. 
The results of wind tunnel tests are also compared with the 
results of some free-flight data for the same wing-body 
combination. 
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CHAPTER I 


INTRODUCTION 


I . 1 Gener al ; 

The problem of wing-body interference is important 
particularly in the case of low aspect ratio v/ing-body 
combinations. The lift and lift curve slopes of low aspect 
ratio wings are generally influenced by the presence of 
bodies. Even in the case of large aspect ratio wings two 
dimensionality of flow on the wing strictly exists only near 
the semi mid-span region. This influence of bodies on wings 
can be expressed in terms of interference coefficients. The 
term interference is used here in the sense of a change in 
fluid dynamic characteristics of a body due to the presence of 
another body. Theoretically the problem has defind a complete, 
and more so, a simple general solution useful in practical 
cases . 

The aerodynamic characteristics of an indivisual 
aomponent of an aircraft are generally known in more certain 
terms. However, the interference of one part on another is 
the most difficult to estimate as they are not easily amenable 
to mathematical analysis, and it cannot be ignored either. 

Many theories and their variations exist for the 
wing alone, the body alone, and their combination 



2 . 


in all the flow regimes. Results of experimental 
investigations for the same theories are seldom available 
to compare. This is more so in the particular cases of 
wings, bodies and their combinations. Evaluation of the 
interference effects experimentally, is the best method 
available to get more concrete results and provide a test 
of the extent of validity of the various theories . 

1.2. Subsonic and Supersonic Theories for Wings, 

Bodies and their Combination, 

The simplest three-dimensional wing theory in 

subsonic flow is the lifting line theory which is applicable 

to unswept wings of high aspect ratio. (AR^6)< The 

solutions of the integro-differential eq.Uation are given 

by Glauert^ , Lotz^ , Miilthopp^ , Gates , Hildebrand and 

Sears^. The lifting line theory is not sufficiently 

accurate for wings of AR-^ 6 or sv/ept wings, if the sweep 

7 

angle > 20°. In such cases Weissinger's theory is more 
appropriate. This^,a first approximation to the lifting 
surface theory. On similar lines lifting surface theories 
due to Blehk^ Ralkner^, H.Schlichting and H.H.H.M. Thomas , 
are also available . In these lifting surface theories 
the Blenk’s theory is the most exact formulation. The 
next in crder of approximations are the theories of Ralkner, 
Schlichting and Thomas. The Weissinger theory cpmes in 
the end. Be Yoig and Barling^ , have given empj-rical 
corrections which bring the seven-point Vsfeissinger method 
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into better agreement with experiment and with the results 

of more accurate calculations based on the 15-point 

12 

Y/e is singer method. Diederich and Zlotnick have given 
sets of aerodynamic influence coefficients based on the 
15-point YYeissinger method for a wide class of wings, they 
have also given the actual spanwise loading distributions on 
the wings and for various power law, twist distribution along 
the span, both symmetric and antisymmetric. In Weissinger’s 
theory it is assumed that the chordwise center-of-lift is 
at the i-chord point. This becomes less and less valid as 
the aspect ratio falls below 2,0 and approaches zero. 

(It is at 20% chord for a rectangular wing of aspect ratio 

1 . 6 ). 

Relatively simple closed form lifting surface 
solution for lift distributions can be obtained for a wide 
variety of V'/ing planforms and angle of attack distributions 
in supersonic three-dimensional flows. This is possible I 
because of the integral equations of linearized lifting 

surface theory may be inverted by a number of special 

15 14 15 

techniques . Harmons and Kainer, J. deal with the cases 

of linearized potential fields and spanwise lift and moment 

distributions for tapered wings with supersonic edges linear, 

quadratic, and cubic twist distributions both symmetric and 

antisymmetric along the span. Bvvard * s ^ ^ method is 

useful in analysing subsonic edges adjacent to a supersonic 

Q- 

edge such as occurs at the tip of,, rectangular wing. 
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This method is not applicable to wings without supersonic 

leading edge as for example delta wings and other highly 

swept back wings. Another planform not amenable to 

Sward's method is the swept wing with subsonic leading 

1 7 

and trailing edges. Works of Cohen gives solution for 
the lift of flat plate swept wings with subsonic leading 
and trailing edges taking into account tip effects. The 
results of this work are often useful in assessing the 
importance of various features of wing geometry on load 
distribution. Numerical methods are available due to 
Etkin and Woodward , and Neuringer These methods 
involve division of the wing surface (and sometimes 
certain regions of wing surface) into grids of 
parallelograms. The contributions of each cell to the 
potential at any point can be summed numerically for 
arbitrary angle of attack distributions over the surface. 

The bodies in fluid flows can be simulated by 

singularities like line sources and line doublets. 

?/heia these singularities are distributed along the body 

axis, they represent the ef feels of body volume, 

incidence and camber. The characteristics of these 

20 

singularities are given in Woodvi?ard's work. 

Bodies at an angle of attack can also be represented 

21 

by circular source rings of varying strength . An 

extensive compilation of theories existing for slender 

22 

bodies are given in the book of J.N. Nielsen 
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There are many subsonic, and supersonic flow 
theories which are available for wing body combinations. 
Subsonic flow theories s 

The mutual interference of a fuselages and wings 

of high aspect ratio in subsonic, incompressible flows have 

23 24 25 

been studied by Lennertz , Wieselsbergh , Pepper and 
Multhopp^^. 

A method for circular bodies involving the use of 

a number of horse-shoe vortices on the wing and their images 

in the body has been suggested by Zlotniok and Robinson ' , 

using conformal mapping to transform the cross-section in a 

28 

TrefXta plane to tnat of a flat plate. Spreiter ana 

29 

V/ard obtained a closed form solution for the case of slender 

wings mounted on cylindrical bodies. The same case with 

30 

elliptical body section has been studied by Nonweiler and 
31 32 

Kahane Mirels has studied the problem of wing-body 
combinations with a gap at the wing-body juncture on similar 
lines . The results of his study are unrealistic as he 

*7 'Z 

neglects viscous and compressibility effects. Lawerence'^ 
gives a theory for low aspect ratio wings with straight 
trailing edges centrally mounted on an infinite ciraolar 
cylinder. This theory holds good even for bodies which 
are not slender for which it was originally formulated. 

The same theory has been extended to arbitrary 

body cross-sections by Lawrence and Plax"^ . (Lagerstrom, 

35 36 3*7 

Graham'^-^ ) jMoskowitz'"^ and (Adams, Sears have studied 
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the lift on after bodies in presence of wing trailing 
vortices. The case of an after-body of decreasing cross 
section area in the presence of wing wake is also studied 
by Adams and Sears. A first approximation for slender 
after-bodies behind non-slender wing-body combinations can 
also be obtained from Ref. 34,36. 

Supersonic Plow Theories; 

Problem of very high aspect ratio wings mounted 
on a conical body considering the interference effect as 
simply the increased lift on wings due to the increased angle 
of attack produced by the upper surface flow around the body 
has been considered by Kirkby and Robinson"^ . Ferrari ’ 
gave an iterative procedure v/hich may be shown to converge 
starting with the body alone and the. wing alone solution and 
procfding to correct the solution in such a way that at each 
stage of the iteration process, the boundary condition' on 
either the wing or the body is satisfied exactly while the 
boundary^ condition on the other is violated. Since the 
process is convergent, a sufficient number of terms will give 
the desired accuracy. Calculations have been done upto 
only the first two steps. His results for the first 
iteration for high aspect ratio rectangular wings mounted 
on a pointed body of revolution show fairly good agreement 
with experiment . 

The problem of rectangular wing and infinite 
cylindrical body of circular cross section at zero angle 



to the main stream on similar lines as those of Ferrai were 

41 

obtained by Nielson . The basis for his formulation of the 

Ap 

problem v/as the work of Lagerstrom and Vandyke . Here the 
basic wing passing through the body gives one part of the 
solution. The flow set up by this wing violates the boundary 
conditions of zero velocity normal to the body. A second part 
to the solution is then obtained which cancels the flow 
normal to the body without at the sametime produc;t5ing any 
new velocities normal to the wing surface,. Thus the 
iteration is terminated in two steps. Nielson solves the 
this boundary value problem betaking the Laplace 
transform of the governing equation. He was able to invert 
the transform of the solution only in terms of a new set 
of functions, one for each terra of the fourier series. Only 
the first four of these functions were tabulated. 

A unified theoretical approach in the field of 

wing-body combination has been tried by many people. Among 

20 

them the work of Woodward's is of considerable importance. 

He deals v/ith the problem by considering the singularities 
simulating the v/ings , bodies and their geometrical 
characteristics. The potential functions for the singularities 
are made to satisfy the Prandtl-Glauert equation along with 
the boundary conditions. As the theory is linear, the 
indivisual potentials can be added and the solution can be 
obtained for the complete combined flow field. To get 
the pressure coefficients for the fields the linearised 
Bernoulli’s equation for steady compressible flow is used 
along with Prandtl-Glauert equation. 
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1,3. Pres ent W ork ; 

In this thesis, the wing alone case is solved by 
Weissinger's method. Aerodynamic characteristics are 
calculated for a class of wings. A rectangular unswept, 
low aspect-ratio, three-dimensional wing with no twist 
and taper, at different angles of attack is studied both 
theoretically and experimentally. Results are compared 
and discussed. A slender body of revolution is studied 
for its aerodynamic characteristics both theoretically 
and experimentally. Aerodynamic characteristics have been 
measured experimentally for a particular cruciform wing- 
body model both in free-f light and in the wind tunnel. 

The results are discussed and compared. 

The instrumentation required for this program 
of studies in wing-body interference has been developed. 
The instrumentation consists of a pair of strain gauge 
balances, a calibration rig for the strain gauge balances, 
models of wing and body of revolution with static pressure 
holes on them. This instrumentation was designed and 
tested out as a part of this vrork to ena.ble a continuing 
program in this field. 



CHAPTER II 


II , 1 . Pprmu.lat_ion_ of^ 

In the following discussion we shall assume 
the wing to be thin £ind the body to be slender, so that 
the perturbations caused by the presence of the wing and 
body combination are small compared to the free stream 
velocity. Plow field due to the wing and the body can be 
represented as equivalent to the surface distribution of 
vortices, over the wing and the surface distribution of 
sources over the body. 


These singularities represent appropriate 
discontinuities in the u,v, or w velocity components 
across the v/ing or body panels and along the body axis. 
The linearized potential function must satisfy the 
Prandtl-Glauert equation. 




The strengths of the singularities are found by 
satisfying the tangential flow boundary conditions at the 
selected control points to get the given configuration. 
The shrouded part of the wing by the body is assumed to 
be of zero camber, zero thickness and zero incidence. 

The theory applies to both subsonic and supersonic flovi/s, 
since the linearized potential function and the above 
Prandtl-Glauert equation hold for these flows. 
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The linearized differential equations for the 
purturbation potential of the wing alone and body alone 
are 



For the combined flow field of a wing and a body 
the total velocity potential is not a simple sum of two 
individual potentials . But the combined potential should 
yield a flow which satisfies the boundary conditions on 
both the wing and body simultaneously. Y/riting down the 
combined velocity potential as 





Many approximate mathematical theories have been 
put forward for the explicit expression of the interference 
velocity potential function.^^’^^’^'^’'^^’^^' Interference 
potential is written as 



.( 5 ) 



where 


(p 


W: 


B.- 


0 


Therefore 


= interference potential due to wing on body 


interference potential due to body on wing 
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(b = -fCp 
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I ^ 
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V/here 
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MW 


1 . ^ 


= the potential function of wing in presence 
of body 


Ha 


tp ' = the potential function of body Inpreaenoe 

of wing 

These three potentials individually satisfy tnt 


equation (i). 


^ ON r 

Q- M ) 4’bcx'^ 4^^= 0 


( 7 ) 


MB hS mb 


( 8 ) 



( 9 ) 



Equation 9 = Equation 7 + Equation 8. 

Since Y is a constant, the second derivative 
of it’s potential is zero, and therefore these potential’s 
may be regarded as the perturbation potential. These 
potential functions are evaluated by appealing to appropriate 
theories mentioned in Chapter I, for the cases of wing alone, 
body alone, and their combinations. The linearized pressure 
coefficients are obtained by linearizing Bernouli’s equation 
for staady, compressible flow, 



vjhere u = (p(bgis "tlie x-perturbation velocity. As the v/hole 
theory is built on small perturbations; it is assumed that 
the same linear solutions also hold good even in the 
neighbourhood of the wing-body junction. 

The pressure discontinuity across the upper and 
lov/er surfaces of the wing panel is given by 



lo\¥er 



upper 
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where hu is the velocity difference on the two sides of the 
wing. The perturbation velocity components are given by 



( 12 ) 


11.2s ^ng_ Alone ; 

(i) Yortex Theory; Y/elssinger 's theory can be considered as 

an extension of lifting line theory or a first approximation 

to the lifting surface theory. The induced angle measured 

at the chord point of the wing is considered to be that 

arising from the bound vortex at the ^--chord line of the wing 

and the trailing vortex and equals the geometrical angle at 

the chord point. The significance of the -f-chord point is 

that in two dimensional flow this is the centroid of the 

function „ / which appears in the ^Arnh’s 

' b i- X 

integral for the lift. In fact for an airfoil section with 
angle of attack plus parabolic camber, the exact two- 
dimensional lift is obtained by considering the entire airfoil 
as a flat plate at an angle of attack equal to the local angle 
of attack of the chord point of the airfoil. 

The integral equation of the Y/eissinger theory may 

46 

be written for swept v/ings as , 




for y>-0; 







^'^±C>Sy>Cp +- 



+ 2+«w 

2 t CLr) (p -^Lc 

>(yi <o) 


( 14 ) 



( 15 ) 
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The coordinate systems -are given in lf'ig.2.1 . 

The boundary condition is satisfied at the f-chord point i.e. 

X = I'VI ta,n(pi- ^ c (16) 

The complete procedure of getting the solution is explained 
in Ref, 46 . 

(ii) Source-Vortex Theory The effe«:ts of wing thickness- 
are given by planar source distributions, The wing camber, 
tv/ist and incidence effects are given by a planar vortex 
system. These singularities can be evaluated in terms of the 
potential function which satisfies the Prandtl-Glauert 
eqn. (1) along with^boundary conditions. The potential 
for the planar singularities can be ^written down as i 

+ VI ( "f. U) 

\/ C i - H G)C y -\}+ zf} 


Cl), 


K 


ri 







k = 0.5 

for 

1 

ll 

0 

for 

M> 1 


( 17 ) 


where 





IL 



represents a constant source distribution 


f 


4 



represents a linearly varying source distribution 


X = Xp- 
^ ^ p ' k 

7. = Zp 


The wing thickness is represented by constant source panels, 
and linear source panels. The wing lift and interference 
is represented by constant pressure vortex panel. 

The integration limits are defined by the panel 
boundaries . For one corner element the potential function 
is given by, (Fig. 2.3 gives the coordinate systems). 


W 

CD = 
k 


k 

2 


A U 
A: 







( 18 ) 
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where 


2AH= pressure discontinuitj across the panel 





la 

^ = y 

cl = 

y' - x -^1 

a' = Lx + Ciy3V 

f ' = y^x- C 

c{' = (iTyy^ ^ 

The potential at a point P(Xp, Yp, Zp) 
is given hy 

W 1 i ^ j ^ 1 ^ 

C^(p)= '4>, - ^ Ty 


where K = 1,2, 3, 4 is the potential due to one 

earner point eleraent of the constant pressure discontinuity. 
The partial derivatives of eqn.(18) gives the perturbation 
velocity component due to each corner element of the 
constant pressure panel. 
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k Aix 
ITT 

- C F3+ F4) 

Vk = 

K AU^ 
2_T 

-^L (F3-f 

Wk = 

2-Tr - y F4 J 

The perturbation velocity components induced at a point 
p by a finite panel may be obtained by super position. 

Up = 


a3 4 d/, 

^p)= 

V. - 

I 


K9(f> 

- 

1a3p_ 


The pressure coefficient Cp is obtained from eqn.(IO) 
and (11). 


II. 3 Bpd^^Alonej 

(iii) Source distributions of varying strength on a circle: 
When the axis of a body of revolution is inclined to the 
main stream, the solution can be obtained by using on the 
surface of the body singularities composed of source rings, 
the strength of which varies round the circumference 


according to, (Pig. 2.2). 
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'j Sl'T) Cp 


where p is ihe a’^imu'thal angle & 

The velocity induced hy an element of such a ring at any 
point (x, T, <p) is given hy 


4- IT Rj 


(21 ) 


The velocity components are obtained by integration 
round the circumference . These are given by , 

sxT- ^c«)- 


rWY ' 




Cp 




2 / + ^ 
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r 2 7>c^) - f w 1 ' 
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( 22 ) 


.....( 23 ) 
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-15 = Qy~,tP) 
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The down-vmsh velocity is given hy 



v/he re 
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( 25 ) 


where the value of K(k) E(k) and D(k) are given hy 

I I ■ I 

K(k) = I ■ clc< 




\y ^ C 2. } 

K O u 'n c\ 


E(k) = 


rTT/i 


\/ 


- olo( 


<0 


T/: 


D(k) 






-d.o^ 


...( 26 ) 


and these are elliptic . integrals. The value of k is 
given by k = / ^ 7 ^ 


where x and r non-dimensionalized with respect to the 
source ring radius r'. (Pig. 2. 2) 


The body is driet-a?iburt©4 byjiurnber of source ring 
along the axis and their induced velocity fields are 
calculated from above equation. The distribution of 
pressure coeffieient over the body can be obtained by, 




2 ,U+ -0 




gfeCA) 
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(iv). The effect of body volume, incidence and camber are 

simulated by line sources and doublets distributed along 

the body axis. The strengths of these singularities are 

determined such that the resulting flow is tangential to 

20 

the surface at each control point. 


The potential which satisfies the eqn. (1 ) for 
the line singularities can be written down as 

^>1 ( cJ'f 


O 




j 


for line sources. (Pig. 2.4) 


& ny 


doublets 


.(27) 


Line source; The potential due to a finite line source 
of unit strength located along the X axis between X = 0 
and X = 1 is obtained by integrating eqn. (27) for n = 1 
bet’ween = 0 and - = . In subsonic flow the integral 

can be integrated directly, but in supersonic flow the 
integral should beCvaluated separately in each of the 
three regions bounded by the Ilach cone from the origin and 
the Mach cone from the tail. (Pig. 2.4). 

In either case the result is given by 



Rc cij — dp - X^/Cog- 


^1-7- d | Nj 

^- 7- <7 2 . J 


. . . . (27a) 


where 


^1 \/ 


d 


2 ■ V + c I - 




The velocity components are given by 
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Line doublet.- For line doublet of unit strength located 
along v/ith X axis is obtained by integrating eqn.(27) for 
n = 2 between limits ^ = 0 and ? = - The potential 

J J 

function valid for both subsonic and supersonic flow is 
given by 
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d.^ , dg are same as before and e is raeasured from the 


^ ~ Y (horizontal plane), axis. The corresponding velocity 
components are 
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II»3. Wlng-bod/_ Ooiabinatipn^ 

The potential function for the combined flow field 
of a wing body combination is obtained by using eqn. (27) 
and (17)* A typical wing body combina.tion is given in 
(Fig, 2. 7). fbe v/ing and the surface of the body are 
subdivided into a large number of quadrilateral panels. 
These panels define the boundaries of planar singularities 
used to simulate effects of wing lift, thickness, and ?/ing 
body interference. Additional line singularities located 
along the body axis are used to siiaulate the lift and 
volume effects of the body. 


The velocity potential function in the present v\fork 
is given by replacing the wing by pressure discontinuities 
across the horizontal panel and the body by a distribution 
of source ring singularities. (Pig. 2. 8). 


The potential function for constant pressure 
singularity is given by eqn.(l7) and the corresponding 
vertical components of perturbation velocity by eqn.(19)» 
Telocity potential for the wing and the corresponding 
induced velocity may be written in the form 
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The down-?/ash velocity at any point P(Xp,Yp,Zp) is given by 

V\/ p ) - \/\l I _ lA/ ^ 4. ^ 

= Cw Li~,w .....(29) 

Vvliere 

+ F^^-UFrFs] 

Ak = 3 

+ [( L> I - M") F^- L(f, - i-'s-) F, 

I 

.(30) 

The graphic-al ■■>■«?. nmabrical integration of the 
velocity field for the body the numerical value of the 
potential function is obtained. When the azis of the body 
is inclined to the free stream a varying circular source and 
distribution is taken as the @iiigul.ftrity . This simulates 
the effect of wing interference on the body. The expressions 




2S 

for the velocity components of a source distribution of 
varying strength on a circle are given by eqn. 22,23,24. 

The resultant, vertical component of velocity due to source 
ring and vortex distribution is given by algebraic sum of two 
components . 

Taking n pressure panel singularities on the wing 
and m source ring singularities on the body we have n 
control points on the wing and m control points on the body. 

The total dovm-wash at any control point on the wing 
is given by algebraic sura of the velocity due to v/i 2 ig 
pressure singularities of eqn.(30) and those due to source 
ring singularities of eqn.(24). 
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(51) 

The down-wash at any of the body control points is given by 
T1 

5._ ^ B 5 ^ mnB ^ Jf' 


( 52 ) 



The down-wash at (m+n) control points is given hy 
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= surface slope at the (m+n) control points 

u 

The system of simultaneous linear equations are 
solved"^*^ by Gauss-Jordon method to satisfy the tangential 
flov/ condition. 


The new potentitil function for the wing in the 
presence of the body is given by equation (34) 
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= pressure discontinuity across the panel 

The value of the potential function at any point 
P(Xp, Yp, Zp) is given by 
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The value of the interference potential due to body on 
wing is given by eqn. (6). 
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Body potential due to source ring is 
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The values for u,v,w, 73,0.. ,'}% are given by 22,23,39, 
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A case of rectangular wing vifith variable sweep 
angles over a body is investigated, for the effect of body 
over spanv/ise lift distributions on the wing. 


CHAPTER III 

EXPERIlffiNTAL JIETHODS USED 

1 1 1 . 1 In_t r od u p t i onj 

The schematic block diagram shown in Pig. 5.1, 
reveals the procedure generally adopted to evaluate the 
coefficients of force, nonents , interference, etc. related 
to an aerodynamic model. The various aerodynamic shapes 
of bodies on which the present studies have been made are 
shown in Pig. 5.2. Design of experiments has been done 
keeping in viev/ the possibility of testing a class of wings, 
a class of bodies, and their combinations by using same 
set of facilities. Most of the facilities have been 
developed specially for these experiments. The table 5A 
shows the details of the facilities and their versatility. 

III. 2. l^esign of Experiments s 

(i) Wind Tunnel? A rectangular low aspect ratio wing with 

no twist and taper was manufactured. The section 
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characteristics are as per NACA 66^012 airfoil . The same 
mould can be used to fabricate wingsof different sweep' 
and aspect ratio shown in Pig. 5.2(b), (c). A slender body 
shown in Pig. 5. 2 (c,d) ?/as fabricated. The body has 4 
sections A,B,C and D) as shovm in Pig.5x2(j) and photo 
no. is By suitably changing the parts A,B, C and D 


a class of bodies can be obtained v/ith differnet nose 
shapes cone angles and slenderness ratios. By suitably 
ch3:.nging parts of Big.J.SCj) and wing sections fig. 5.2^0 
and(c) we get v/ing-body combinations as shown in fig. 3. 2 
and (h). By manufacturing different v/ings and 
sections A,B,C,D of Pig.3.2(j) we can generate a large 
number of wing-body coiabinations of practical interest 
The size and strengths of the models depends on the speed 
range and working sections size available in the wind 
tunnels. A strain gauge balance which can measure 
5 components was designed out of S97 steel (B.S.S.) for 1 
fairly heavy models. Another 6 ^ components strain gauge 
balance was designed for light models out of L83 light 
alloy , (B. S .3 . ) keeping in view the requirements of free 
flight testing and the existing lov/ speed tunnel work. 

A model mounting system v/as designed to take up loads 
on the model by using 5 component balance mounted 
external to the tunnel, fig. 3. 2 (k)(s) and (t) shows ■ 
various arrangements. A calibration rig was designed 
to calibrate both the 5 component strain gauge balance. 
The capacity of the rig was kept fairly high to cater to 
the needs for calibrating high load range strain gauge 
balances . (for example balances suitable for transonic, 
shock loads, supersonic and high subsonic ' flows ,) . 
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The coordinate systems, for which this calibration rig 
is suitable for calibration of strain gauge balances are 
shown in Pig. 3. 3. 

¥/ith these facilities it is clear that a class 
of wings , bodies , and their combinations can be studied 
at low speeds . 

(ii) Pree-Plights The model shown in Pig.3.2(e) was 
mounted externally on a five component balance. The whole 
assembly v/as mounted on the "nose boom" of a glider also 
termed the" aerodynamic platform. Pig. 3.4 shows the 
sketch of the arrangement. The Rohini glider v/as towed by 
an powered aircraft 'Piper’ . The level of flight • 
for glider v/a.s kept well above compared to the airplane 
during tov/. When sufficient height was reached the glider 
was released and v/as made to sink at a desired constant 
speed. Then the force measurements v/ere done on the model 
for various angleSof attacks, (adjusted from the cockpit). 
The read-out system was a simple portable budd strain 
indicator. The results are discussed in Chapter V. 


CHAPTER IV 


PROCEDURE 


IV.1, Intro duction; 


The whole set of experiments was carried out 
in the (2'x3'x5') 3-D test section (Photo 1) of the 
lowspeed wind tunnel of the Department of Aeronautical 
Engineering. As part of the test requirement, information 
regarding the flow characteristics in the test section was 
obtained through a detailed calibration of the tunnel. The 
strain gauge balances designed to measure the force systems 
on the models, were calibrated. Experiments on a rectangular 
wing (Photo 3»7), a body (Eig.3.2 (d) , (t ) , Phoiio Hos,2,17) 
and a cruciform wing body combinations (Fig. 3. 2(e), Photo 
No. 1), were carried out. The cruciform wing body combination 
was tested both in the wind tunnel and in free-f light. 
Qualitative studies are made of the trailing vortex system 
for the rectangular wing model. (Phono No. 8, 8). 

IV. 2 Tunnel Calibratio n ; 

Calibration of the flov/ in the tunnel test section 
is carried out for the following objectives ; 
a). Velocity distribution in the jet at various 

cross sections . (Pig. 4.1 ) . 


b) . i’low inclination of the jet. 

c) . Turbulence level. i < 3 , 2^2.'’ ? (h.) j 

d) . longitudinal static pressure gradient. 

a). Dynamio head survey at various cross section 

were carried out by using a pitot-static tube (Fig. 4 . 2). 

The dynamic pressure distributions in the rectangular 

test section are shown in Fig. (4-3 A,B,C,D). The region near 

the walls was avoided as the probe was subjected tot 

strong boundary layer effects. The laminar boundary layer 

grov/th along the test section wall was calculated by using 
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the expression for a flat plate. 



Value of _ was equal to 0.01213 inch- For a distance of 
66" length. It was qualitatively observed that it was less 
than 0.2". (by using a thin total head tube and a stethoscope) 

B). Flow inclination in the jet was measured by 
using three types of yaw probes(Fhoto ho.10). All of them 
gave very nearly the same result. The flow inclination in 
the tunnel is 0,3 degrees, tovmrds bottom wall of the test 


section. 


m 



G). The turbulence level in the wind tunnel section 
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was calculated by using a sphere The turbulence factor 

in the tunnel is 1.35. 

D). The observed static pressure gradient along the 
test section is plotted in |'ig.4-.4. 


IV. 3 . Str ain Gau ge Balance_ hesigns : 


It is generally necessary to have a number of 
balances for different load ranges, and of different sizes 
for the same load range, for testing various type of models 
in a wind tunnel. There are many types of wind tunnel 
balances available in concept for design and fabrication. 

Each of them has it's ovm advantages and special characteristics. 
Survey of the litrature for concept revealed the characteristics 
which are tabulated in Table IV. A. The resistance type of 
strain-gauge balances are in common use as they are easy to 
design, fabricate, and simple to operate. The two strain 
gauge balances are designed in conformity v/ith the following 
salient oonoeptsS^. 55 , 56, 57 ,58,59 . 


1 ) . The material chosen has the following 
characteristics i 

a) , good elastic properties 

b) . easy machinebility 

c) . minimum y/arping under heat treatment 
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d ) . compatibility with the strain gauges in respect 

of .expansion coefficient and adhesion qualities 

e) . good strength and fatique properties 

f) . lov/ thermo-elastic coefficients 

2) . Deflections under load to be minimum to 

give 

a) , lov^ interactions 

b) . small clearance required between model & balance 

c) . elimination or reduction of model attitude control 

d) . high natural frequency 

3) . The measured strains were kept within 0.003 
to 0.001 inch/inch to keep the temperature effects to a 
minimum. 

4) . Low stressed areas under gauge grid and low 
stress area under lead-wires. 

5) . Gauging area accessible for ease of locating, 
mounting, and servicing of gauges, but gauges should be 
protected to prevent damage from fouling under load or 
handling. 

6) . Precautions in choosing machining tolerences 
to make sure the interaction due to dimensional differences 
or angular misalignments be kept to a minimum. 
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7 ) . Strain gauges of the same bridge are exposed 
to equal temperature . 

8) , Strain gauges locations to be placed whese 
unwanted asymmetric high stresses do not exist. 

9) . The spring constants are the same for both 
positive and negative loads to eliminate two sets of 
calibration data. 


The drawings of the balances giving physical 
dimensions are shown in Pig. 4. 5 and 4.6. The specification 
and load ranges for the balances are as follov/s^ 


a). Piye-Corappnent_^a_lance_; (Photo No.11 ) 

Material Steel S.97 (B.S.S.) 

Physical dimensions 2"x2" square section of 9" length 

Weight 1.2 kgs. 


Component full load ranges 


Normal force 
Pitching moment 
Side force 
Yawing moment 
Chord force 


100 lbs. 
100" lbs. 
1 00 lbs . 
100" lbs. 
15 lbs. 


Cp travel v/ith the model allov/ed from the centre of the 
balance is + 1 " . 



^0 


3 ) . S ix-c o_mp one n t Balance ; 

Material 

Phy s i c al a ine ns i ons 

V7eight 

Component full loEd ranges 
Normal force 
Pitching moment 
Side force 
Yawing moment 
Chord force 
Rolling moment 


(Photo No.12,13) 

Light aircraft aluminium alloy 
L.83 (B.S.S.) 

0.75" diameter cross section 
Length; 7-5" 

C.75 lbs. 

50 lbs. 

50" lbs. 

50 lbs. 

50" lbs . 

5 lbs. 

50" lbs . 


Cp travel with the model allowed from the centre of the 
balance is + 0.75". 


IV.. 4 Calibratic 


This equipment helps in loading the strain gauge 
balances and calibrate them, (Phcto Nos .1 4, 1 5 , 1 6, 1 9 ) dhows 
a' details of the equipment. The balance along with it’s 
calibration fixtures, (Pig.4.7,I'Boto Nos .18,1 5,16) are fixed 
on the rig and taken through the axis systeS^shown in 
Pig. 4. 8. At each position of the coordinate system shov/n 
we will be loading normal force pitching moment , side force , 
yawing moment, rolling moment and drag respectively. 
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The loading v/eights were purchased commercially and calibrated to 

an accuracy of + loading is done through cables. 

The sting end of the balances are to be of one standard 
dimensions, so that they can be quickly assembled on the 
Eig. (Pig. 4.7). Othervi/ise adapters have to be manufactured 
before starting of work. The loading range capability of the 
rig is given below: 


Normal force 
Pitching moment 
Side force 
Yawing moment 
Rolling moment 
Drag 


1000 lbs. 

1 000" l^bs -. 
1000 lbs. 

1 000" !-^bs . 
500 fjjlbs . 
100 lbs. 


Though the present balances do not use the full 
capabilities of the rig, it is built keeping in vj.ew future 
possible needs as already explained in Chapter III. 


A static load deflection survey of a balance can 
be conducted by using the rig. This information will be useful 
in making corrections for angle of incidence v\fith models in 
the airflow. 


I V . 5 . C alibrati on of v /in d tun nel B alance; 

This phase of the work is very important in making 
the effort of many man days of design, fabrication worthwhile 
in terms of dependable aerodynamic measurements. Calibration 


4,2 


of the balance has the following object ivesf^ 

a) , to proof load the balance 

b) . to determine calibration slopes for each 

component 

c) . to determine component sensitivity 

d) . to determine balance interactions 

e) . to determine deflections under load 

f) . to check repeatability of load data 

Although in the balance design the greatest possible 
emphasis is placed on eliminating sources of non-repeatable 
or indeterminate errors, large interactions may appear, due 
to dimensional limitations and manufacturing inaccuracies 
on balance design. However, these interactions are 
repeatable and can be accurately determined, so that the 
balance output data can be corrected with acceptable accuracy, 

To load the balance accurately a fixture which 
simulates the model to be tested is fabricated . (Photo No. 15, 16, 18 
Pig. No. 4. 7). This v/ill transmit loads to the balance in 
the same manner as the model during a wind tunnel test . This 
fixture should be attached precisely as the model is attached 
to the balance, and should have reference surfaces which will 
remain fixed with respect to the balance reference, as the 
balance, being loaded. ^ 
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The output from an indivisual bridge is due to it 
own component loading plus that due to loading of other 
components. The first one is called as "component output", 
and the second, as "interaction output". Thus for example 
the normal force output can be expressed as 0^ = ^ 

is the normal force output. This can be written as 


© 


1 


where 


I ^7]) 


|\! = normal force component output 

X| (^) = normal force interaction output 


We assume here only linear interactions. Writing down in 
general terms the output of each indivisual bridge as below 
will simplify the procedure for evaluating calibration 
constants for a balance. 


Si ^ -f- CV2 2X 4- 4 

SS + 4 +-4^^-^ f 03^4) 

SS + !JL^,x 0^45 Y t ■)- 44^^ 

Z Z %_2X ->■ Q 4 Q^yfri 4 

■+ CI 43 Y •+ ^4^ 


This can he written dov/n further as 


L©J = Ca] l f] 




<^22 ^23 

^3) ^32 ^33%k%Ac 

«« 

^1 «a <^3 c^4. <?Cs-<^ 



= output Matrix 



coefficient Matrix 



Force Matrix 


Each diagonal tern of the Matrix 



is the 


component output slope. The other coefficients of the 


Matrix {' f are the slopes due to interactions 

The aim of a good balance design is to get (i) the 
coefficients other than the principal diagnjnal coefficients 
equal to zero in the Matrix jand (ii) to get as high a 

value as possible for the principal diagonal coefficients. 




This means , a sensitive, interaction free balance. The 
criteria is then of diogonalizing the Matrix i/'j 
The terms Q--! ^ > i / D should be expressed in terms of 
the geometry and tolerances of manufacture of the balance 
indimensionless form. In approach in this direction may 
contribute to the design of balances . 


The force Matrix in terms of inverse of 
coefficient Matrix and the output Ifetrix is 


[i~J= [Aj j_ e”| 
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The Matrix ^ fixed property of the balance, so that 

for any measured strain Ivlatrix the corresponding force 
Matrix is obtainable. The results of calibration of the 
five component balance is plotted in Pig. 4. 3- fbe six 
component balance is fabricated and strain gauged. As 
the wiring is not yet ready the calibration data is not 
available. All the experiments were carried out by using 
five component balance. 


IV, 6 . Wi ng IJode l; 

A rectangular symmetrical wing of aspect ratio 

2,5 with no twist and taper vms fabricated. (Photo IJo.3j7)* 
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The airfoil section is of NACA 66^-012 data . As the wing 
has a symmetrical cross section the top half surface to 


the left and hottora half surface to the right is tapped 
for pressures. A total of 12 chord positions pressures 
can be tapped one at a time. Each chord contains 15 
pressure taps . At various angles of attack in the flow 
the pressures were measured and simultaniously, the 
balance output were taken and the results are shown in 
Eig. 4. iO, Th§6odel v/as exited by a step function without 
wind and v/ith v^ind. The time for decay to half amplitude 

iti;" KC 

in both cases . ' calculated and compared. This gives an 

idea of the aerodynamic damping on the wing. (Appendix A). 

I V . 7 • Body Zlo del ; 

A slender body is fabricated in four separable 
parts. (Photo No.l?). The reasons for doing so has already 
been explained in Chapter III. The distribution of pressure 
coefficients along the model is measured at various angles 
of attack in the flow. The same is plotted in Pig.4.V|. 
Simultaneously balance outputs v/ere recorded and plotted. 
(Pig.4.11-). 

Dynamic measurements are done on the body, just 
on the same lines as for the v/ing. The results are in 
Appendix A, The body mounted in the tunnel test section 
is shown in Photograph No. 2. 
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IV. 6. ^iicifoTm_ }lin^-Body^ Pprnb in_at i onj_ 

A slender body v/ith cruciform wings is shown in 
Photograph No.1. This model was tested for force measurements 
using the five component strain gauge balance. First it was 
airborne on an "Aerodynamic Platform" as explained in Chapter III.. 
Then the same model was tested in a wind tunnel. The results 
are shovm in Pig. Zp' 13. 

Dynamic measurements are carried out on this model, 
as in the case of wing and body mentioned above . The results 
are shown in Appendix A. 

IV, 9. Vortex Studies ; 

The trailing vortex system for the refttangular wing 
at negative stalling angle of attack is visualized by using 
a tuft grid. (Photo 110,8,9). The g^id is moved, behind the ?/ing 
and observed trailing vertices. The results are shown in 
Pig. 4. :|.^The followring points were observed in these studies? 

a) . The tv/o trailing tip vortices are symmetrically placed. 

b) . The starting of tip vortex takes place just about 

of the inboard. 

c) . The tip vortex is circular in shape. 

d) . The size of the vortex increases v/ith the distance 

e) . The vortex system moves more and more inboard with 

distance. 

f) . Appears to have an overlapping shedding frequency 

over and above the vortex system 


CHAPTER Y 


DISGPS SIQH Aro COIC LU SIORS 
V.1. I ntroduction; 

In Chapter I, the various theories for wings, bodies, 
and their combination in 3-diniensionaI flows were reviewed. In 
Chapter II, the formulation of the problem was done. 

In Chapter III, the experimental technique used in the present 
investigation is explained. In Chapter IV, the procedures 
followed to obtain the experimental results are discussed. 

The results obtained from such procedures are discussed in 
the following articles . 

V. 2 . General ; 

The flow characteristics in the tunnel test section 
are Irery suitable for experiments of pressure plots, and force 
measurements. The corrections regarding static pressure gradient, 
flow inclination are negligible. The high value obtained for 
the turbulence factor is on hi«gher side as the turbulence sphere 
used for measurement was not either spherical or smooth 
precisely. It can be expected to be much less than 1.35. 

It is necessary to make a more accurate measurement of the 
turbulence factor using a better sphere than was available 
for the present experiments. As for the present investigations 
it is hoped this will not effect the results seriously. 

The strain gauge balance calibration was conducted 
only on the five component balance. The six component 
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balance is now ready for wiring. After wiring the balance can 
be calibrated and used for further investigations . ffrom the 
calibration results of the five component strain gauge balance 
it is seen that the balance measurement are well within 0.1^ 
of the full load for all the arms of the balance. 

V.5. W ing Model; 

It may be seen from Fig. 4. 9 that the chordwise 
pressure distribution agrees well with the theoretical values. 
The descripency in the experimental points on chord stations 5»6 
are due to the influence field of the support system. 

(Photo No. 3 , 1 ) • The spanwise lift distribution calculated by 
Weissinger's method and Woodward's method agree very well with 
the experimental results of pressure plotting. The force 
coefficients obtained by pressure plots agree well with the 
force measurements done by the five component strain gauge 
balance . The only Slscrupen^y occures in the drag coefficient 
plot. The drag coefficients obtained by pressure plots do not 
agree with the measured drag force^one of the reasons being 
that the pressure plots do not include^ the induced drag. 

The trailing vortex system vms studied qualitatively. The 
location and field of the trailing vortices are plotted from 
a survey of the tufts grid behind the wing. A closer 

measurement of the down-wash field of the wing and wing— body 

the: 

combination will be necessary to estimate^ influence field due 
to the interference. 


So 


The dynamic measurements conducted un the v/ing showed 
that it is a highly damped system in the air-flov\f for pitching 
disturbances. From n comparison of the measurements of the 
damping in pitch of the wing alone v^ith and without wind 
(Appendix A). The time for decay to half the reference 
amplitude for the Vi/ind-on case was I/ 3 , that of the wind off 
case. The inertia system of the balance and support are tafceated 
as constant in this case . 


V.4. Bod y Model; 

A body model for pressure plots was tested in the 
wind tunnel. The results shov/ ei good agreement with the 
theoretically calcula.ted values for a varying source ring 
distribution. At zero degree incidence the flov; is axisymmetric 
over the body. This checks indirectly the flow characteristics 
of the test section also. The discrepancy in the hiegher 
angles of attack can be attributed to 

a) . the viscous effects 

b) . the interference field due to support system 

c) . influence of base pressure upstream 

The more predominant of them is the interference of 
model support system. The other contributions can be assessed 
only after further investigations. 

The flow over the cone is smooth at lower angle of 
attack. At hiegher angles the flow may have separa,ted and 
possibly reattached to the body. The conclusions anyway are 
open for further investigations. From an integration of the 


pressures over the body it is seen that most of the lift of 
thE body arised from the conical portion. 


The values of the force coefficients obtained 
for this model from strain gauge balance are also plotted. 

The damping coefficient in pitch is small for this body but 
is found to be ^ that of the wind off case. 

V.5. Cru ciform Wing-body Model; 

V/ith a view to develop a technique of testing 
aerodynamic models airborne, particularly for V/STOl models 
an attempt was made to test this wing -body model both in 
"free-flight " and in the wind tunnel. The same five component 
strain gauge balance v/as used in both experiments, A plot 
of the results reveal that a good future exists there for 
testing airborne models. The scatter in free-flight data 
can be minimised by designing a low interference model support I 

system. This technique eliminates the strong wall interference ; 

I 

effects present in conventional wind tunnels on v/STOl models, ' 

V.5. The wing-b ody C ombination; j 

I 

Theoretical calculations using the method of ‘ 

20 ^ * 
?/oodward's were made for the present wing-body combination, | 

i 

i 

The interference field of the body on the spanwise lift i 

distribution of the wing for sweep angles 0°,45°,60° were 

i: 

calculated. Measurements for comparing these values with | 

1 ^ 

experimental values are to be continued further. ; 

) Clnd -$"'1, , I 

I i. !.T Ai.'-'.'ilR, I 

i Li---' I I 

t d I 

I i: 

_ _ _ _ _ ^ 


V.7» Conclusions ; 


The present work has formed a basis for further 
work in the field. It is possible to investigate many bodies, 
wings and their combinations and generate valuable 
aerodynamic data for wing-body combination at low speed. 
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APmroiji _ A* 

The equotion of motion for a single degree of 
freedom pitching system with no air loa-ding or air damping 
hut v/ith a mechc.nical spring that tends to return the model 
to an equilibrium position is as follov^s; 


.1_0 p -f i . y f -f- K P 


'-f 
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(1)A 


v/ith airs 


0 + pc, + c^-,) e-+(k,+ !<de = 0 (2)A 


Solutions of eqns. (1 )A and (2 )a are 
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Kp- -L 4T 
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pi [O' 6951^ 




C'kr m 






(10)A 


where subscript f is for case of v/ithout air 



inertia of the model about the axis of rotation 
s lugs/ft ^ 

angle of attack of the model to the equiULbrium 
position (usually zero, radians) 
a system damping constant ft lbs/( radians/sec) 
a system spring constant (ft Ibs/radians) 


c = 

k = 

2^ 



aerodynamic "damping moment" per unit value 
of ft lbs/(radians/sec ) 

.. -.or'i ' rtt 

aerodynamic "restoring laeemtn" per unit value 
of ft Ibs/radians 

model pitching frequency cycles/sec. 
period of one cycle of oscillation sec 
time required for the amplitude of model 
oscilla^tion to decay from reference value to 


a value equal to half the reference value sec. 

Ref. F Ki A 1_ 


TABLE A(1 ) 


• ♦ *' 


Prom 

Measurements (Graphs) 

— • (f 


SI*. 

Bo. 

Case iff 1 

■Cvl/*.c.j 


■k 

.^./CyCi 


1 . 

*Wing-body conbinatior 

0.5 

3.34 

1 .8 

0.2 

5.0 

• ♦ « 

2.1 


(Cent are Rocket Model) 







C<=: 0 ° 







2. 

^Rectangular wing 








alone 

0.16 

6.67 

2.24 

0.140 

7.15 

0.16 









3 . 

*Body alone 








(Pressure Model) 

0.5 

3.34 

1 .16 - 

0.35 

2.58 

2.25 


o<- o'" 







4. 





*500 

3.334 

2.00 







4f “* • ”“ i 



* All the models were exited by a step forcing function while 
they were mounted on the strain gauge balance. ^42,^. 
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